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We introduce a class of multifractal processes, referred to as Multifractal Random Walks (MRWs). 
To our knowledge, it is the first multifractal processes with continuous dilation invariance properties 
and stationary increments. MRWs are very attractive alternative processes to classical cascade-like 
multifractal models since they do not involve any particular scale ratio. The MRWs are indexed 
by few parameters that are shown to control in a very direct way the multifractal spectrum and 
the correlation structure of the increments. We briefly explain how, in the same way, one can build 
stationary multifractal processes or positive random measures. 

PACS numbers: 05.45.Df Fractals, 47.53.+n Fractals in Fluid dynamics, 02.50.Ey Probability theory, stochas- 
tic processes, and statistics, 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion 
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Multifractal models have been used to account for scale 
invariance properties of various objects in very diffcr- 
ent domains ranging from the energy dissipation or the 
velocity field in turbulent flows to financial data. The 
scale invariance properties of a deterministic fractal func- 
tion f(t) are generally characterized by the exponents ţ q 
which govern the power law scaling of the absolute mo- 
ments of its fluctuations, i.e., 



m(q,l) = K q l c \ 



(1) 



where, for instance, one can choose m(q,l) — ^2 t \ f(t + 
l) — f{t)\ q . When the exponents ţ q are linear in q, a 
single scaling exponent H is involved. One has Cq = qH 
and f(t) is said to be monofractal. If the function ( q 
is no longer linear in q, f(t) is said to be multifractal. 
In the case of a stochastic process X(t) with stationary 
increments, these definitions are naturally extended using 

m(q, l)=E (\SiX(t)\") = E (\X(t + l) - X(t)\i) , (2) 

where E stands for the expectation. Some very popular 
monofractal stochastic processes are the so-called self- 
similar processes 0]. They are defined as processes X(t) 
which have stationary increments and which verify (in 
law) 



S X iX(t) = X H S l X{t), VZ,A>0. 



(3) 



Widely used examples of such processes are fractional 
Brownian motions (fBm) and Levy walks. One reason 
for their success is that, as it is generally the case in 
experimental time-series, they do not involve any par- 
ticular scale ratio (i.e., there is no constraint on l or A 
in Eq. (||)). In the same spirit, one can try to build 
multifractal processes which do not involve any particu- 
lar scale ratio. A common approach originally proposcd 
by several authors in the field of fully developed turbu- 
lence [| [|, [|, 0], has been to describe such processes 



in terms of differential equations, in the scale domain, 
describing the cascading process that rules how the fluc- 
tuations evolves when going from coarse to fine scales. 
One can state that the fluctuations at scales / and XI 
(A < 1) are related (for fixed t) through the infmitesimal 
(A = 1 — r] with rj « 1) cascading rule 



5 xl X{t) = W x SiX(t) 



(4) 



where W\ is a stochastic variable which depends only on 
A. Let us note that this latter equation can be simply 
seen as a generalization of Eq. (||) with H being stochas- 
tic. Since Eq. (|j) can be iterated, it implicitely imposes 
the random variable W\ to have a log infinitely divisi- 
ble law. However, according to our knowledge, nobody 
has succeeded in building effectively such processes yet, 
mainly because of the peculiar constraints in the time- 
scale half-plane. The integral equation corresponding to 
this infmitesimal approach has been proposed by Cas- 
taing et al. ^ . It relates the probability density function 
(pdf) Pi(SX) of SiX to the pdf G X i,i of In W x : 



Pm(6X) - / G m {u)e- u P l {e- u 8X)du. 



(5) 



The self-similarity kernel Gaz,z satisfies the same iterative 
rule as W\ which implies that its Fourier transform is of 
the form G\i^{k) — G x (k). Thus one can easily show 
that the q order absolute moments at scale / scales like 



m(q, l) = Gi, L (-iq)m(q,L) = m(q,L) 



F(-iq) 



(6) 



where F = In G refers to the cumulant generating func- 
tion of hiW ^,0]. Thus, identifying this latter equation 
with Eq. (0), one finds £((?) = F(-iq). 

In the case of self-similar processes of exponent H , one 
easily gets that the kernel is a dirac function G^l{u) — 
S(u — Hln(l/L)) and consequently ţ q — qH. The sim- 
plest non-linear (i.e., multifractal) case is the so-called 



2 



log-normal model that corresponds to a parabolic ( q and 
a Gaussian kernel. Let us note that if £ 9 is non linear, 
a simple concavity argument shows that Eq. (jî|) cannot 
hold for all l in ]0, +oo[. Multiplicative cascading pro- 
ccsscs § |, [fi], y, [TJ consist in writing Eq. (g) starting 
from some "coarse" scale £ = L and then iterating it 
towards finer scales using an arbitrary fixed scale ratio 
(e.g., A = 1/2). Such processes can be contructed rigor- 
ously using, for instance, orthonormal wavelet bases [ O . 
However, these processes have fundamental drawbacks: 
they do not lead to stationary increments and they do 
not have continuous dilation invariance properties. In- 
deed, they involve a particular arbitrary scale ratio, i.e., 
Eq (|î|) holds only for the discrete scales l n — X n L. 

The goal of this paper is to build a multifractal process 
X(t), referred to as a Multifractal Random Walk (MRW), 
with stationary increments and such that Eq. (|l|) holds 
for alU < L. We first build a discretized version X At (t = 
kAt) of this process. Let us note that the theoretical issue 
whether the limit process X{t) = limAt^o X At (t) is well 
defined will be addressed in a fortheoming paper. In this 
paper, we explain how it is built and prove that different 
quantities (q order moments, increment correlation,...) 
converge, when At — > 0. 

Writing Eq. (||) at the smallest scale suggests that 
a good candidate might be such that 5 At X At (kAt) = 
£At[k] W A t [k] where t At is a Gaussian variable and 
WAt[fc] = e" A '[ fe l is a log normal variable, i.e., 



t/At 

X A t(t) = S AtX At (t) 

k=l 



t/At 

E 

k=l 



£At[k)e' 



(7) 



with X At (Q) = and t = kAt. Moreover, we will choose 
t A t and LOAt to be decorrelated and e At to be a white 
noise of variance a 1 At. Obviously, we need to correlate 
the WA t [fc]'s otherwise X At would simply converge to- 
wards a Brownian motion. Sincc, in the case of cascade- 
like processes, it has been shown JÎ2], [Î3| Q that the 
covariance of the logarithm of the increments decreases 
logarithmically, it seems natural to choose 



Cov{LOAt[ki], WA<[fe]) = A 2 lnpAt[[&i - k 2 



with 



PAt[K\ 



1 



-l)At 



for |fc| < L/ At 
otherwise 



(8) 



(9) 



i.e., the w^t are correlated up to a distance of L and their 
variance A 2 ln(L/At) goes to +00 when At goes to 0. For 
the variance of X At to converge, a quick computation 
shows that we need to choose 

E [w At [k]) = -rVar [ui At [k]) = -r\ 2 ln(L/At), (10) 

with r = 1 (this value will be changed later) for which 
we find Var{X{t)) = a 2 t. 



Let us compute the moments of the increments of the 
MRW X(t). Using the defmition of X At (t) one gets 

E(x At ( tl )...x At (t m )) = El 1 /^ -Et= A i* 

E{e At [k 1 ]...e At [k m ])E (e«A»[ki]+...-u»A*[fcm]) . 

Since e At is a mean Gaussian process, this expression 
is if m is odd. Let m = 2p. Since the est[k]'s are 
(5-corrclatcd Gaussian variables, one shows that the pre- 
vious expression reduces to 



T 2p 



(*S(l) At S(2))/ At < t S(2p-l) At S(2p))/ At 



— y t 

ses 2p fci=i 



E (>£?=i WA *fel) At p , 



where aAb refers to the minimum of a and b and S%p to the 
set of the permutations on {1, 2p}. On the other hand, 

we have E (^e 2 ^=^ At[k ^ = Ui<j P[ k i ~ ■ Then, 

when At — > 0, the general expression of the moments is 



E(X(t 1 )...X(t 2p )) = f^j Eses 2p Io 



*S(1) ^*5(2) 



du^ 



where p(t) = lmiAt^o p A t[t/Ai\. In the special case t\ = 
ti = ... = t 2p — l, a simple scaling argument leads to the 
continuous dilation invariance property 



m(2p, l) = K- 



2p 



p-2p(p-l)X 2 



VI < L, 



(12) 



where we have denoted the prefactor 

K 2p = L p a 2p (2p-l)\\ f dui... [ du p T\\ui-Uj\ 



By analytical continuation, we thus obtain the following 
C g spectrum 



( q = (q - q(q - 2)A 2 )/2. 



(13) 



We have illustrated this scaling behavior in fig. 1. Thus, 
the MRW X(t) is a multifractal process with stationary 
increments and with continuous dilation invariance prop- 
erties up to the scale L. Let us note that above this 
scale (/ >> L), one gets from Eq. ( |î"î| ) that ţ q = q/2, 
i.e., the process scales like a simple Brownian motion, as 
if oj was not correlated, though, of course, X(t) is not 
Gaussian. Indeed, K 2p is nothing but the moment of or- 
der 2p of the random variable X(L) and is infinite for 
large p's (depending on A). Actually, one can show that 
C2p < K 2p = +00. Consequently, the pdf of X(L) 
has fat tails. As illustrated in fig. 2, Eq. (|J) accounts 
very well for the evolution of the pdf of the increments. 
One shows that the smaller the scale /, the fatter the tails 
of the pdf of 5iX{t). 
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FIG. 1: (a) Plot of two realizations of 2 17 samples of two 
MRWs with A 2 = 0.03, L = 2048 and where e A t is (top plot) 
a white noise or (bottom plot) a fGn (Eq. (LLa)) with H — 2/3 
. (b) Log-log plots of m(q, l) of the MRW plotted in (a) (top 
plot) versus l for q — 1, 2, 3, 4, 5. (c) (o) (resp. (+)) : C, q spec- 
trum estimation of the MRW plotted at the top (resp. bot- 
tom) in (a) . These estimations (obtained using the WTMM 
method jla]) are in perfect agreement with the theoretical 
predictions ( ) given by Eq. (|Î3| ) (resp. Eq. (|l9|)). 



Let us study the correlation structure of the increments 
of X{t). Since ( 2 — 1, one can prove that they are decor- 
related (though not independant). Let 



C 2p (1,t) =< \S T X(l)\ 2p \5 T X(0) 



2p 



>, 



(14) 



with t < l. Using the same kind of arguments as above, 
one can show that 



C 2p (l, r) = (a 2p (2p - l)!!) 2 l! +T du,... J/ +r du 



r-l+T 



l ""*P 
4A 2 



$a du v+i- $a du 2pTli<i< 3 <2 P p( u i - u jY X "- (15) 
A straightforward argument then shows that 

2p ((l + t) /L)^ 2 p 2 ~ 2p[ ' ' ~ 2p ((Z - t) /L)^ 2 p 2 ' 



ti 




FIG. 2: (a) Standardized estimated pdf's of hi SiX(t) for 
« = 4,32,256, 2048 and 4096 (from top to bottom). These 
estimations have been made on 500 realizations of 2 17 samples 
each of a MRW with A 2 = 0.06 and L = 2048. Plots have been 

arbitrarily shifted for illustration purpose. ( ) theoretical 

prediction from the estimated pdf at the largest scale (l = 
2048) using the Castaing's equation (^|). 



and consequently for r << l fixed, using analytical con- 
tinuation one expects C q (l,T) to scale like 



C q (l,r)~Ki 



(i) 



2C, 



—lambda q 



(16) 



This behavior is illustrated in fig. 3. 

^From the behavior of C q when q — > 0, we can obtain 
using Eq. (^6|) that the covariance of the logarithm of 
the increments at scale r and lag l behaves (for r << l) 
like 



C (ln) (Z, 



(17) 



Thus, this correlation reflects the correlation of the u>At 
process and is the same as observed in Refs Jî^, [b|] 
for the cascade models. This behavior is checked in fig. 
4. 

Finally, let us note that, one can built MRWs with 
correlated increments by just replacing the white noise 
£At by a fractional Gaussian noise (fGn) 



4^[fc] =B H ((k+l)At)-B H (kAt), 



(18) 



where Bn(t) is a fBm with the scaling exponent H and 
of variance a 2 t 2H , and choosing r = 1/2 in Eq. (|î"c|). 
Then, one can show (after tedious but straightforward 
computations) that the spectrum of the MRW X^ H '(t) 
is 



Cf > =qH- q(q - l)A a /2, 



(19) 



(Q = at scales >> L) and consequently the MRW 
has correlated increments. Such a construction is illus- 
trated in fig. 1 with H = 2/3. Since H > 1/2 it leads to 
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FIG. 3: (a) Log-log plots of C q (l,r) versus l for q = 1,2,3. 
(b) Estimation (o) of the power law exponent C q (l,r) ~ l Vq . 
It is in perfect agreement with the theoretical prediction (Eq. 
©) „, = -A 2 q 2 ( ). 
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FIG. 4: Estimation (o) of the correlation C (ln) (/,r) of the 
logarithm of the increments. It is in perfect agreement with 
the analytical expression (Eq. (^7D) — A 2 ln(-^) ( ). 



a process which is more regular than the one previously 
built. 

To summarize, we have built the MRWs, a class of mul- 
tifractal processes, with stationary increments and con- 
tinuous dilation invariance. ^From a theoretical point 
of view, MRW can be seen as the simplest model that 
contains the main ingredients for multifractality, namely, 
logarithmic decaying correlation of the logarithms of am- 
plitude fluctuations. Moreover, they involve very few pa- 
rameters, mainly, the correlation length L, the intermit- 



tency parameter A 2 , the variance a 2 and the roughness 
exponent H . They all can be easily estimated using the 
( q spectrum and the increment correlation structure. We 
do believe that they should be very helpful in all the 
fields where multiscaling is observed. MRWs have al- 
ready been proved successful for modclling financial data 
p6[ . In this framework, we have shown that one can eas- 
ily build multivariate MRWs. Actually, the construction 
of MRWs, as presented in this paper, can be used as a 
general framework in which one can easily build other 
classes of processes in order to match some specific ex- 
perimental situations. For instance, a stationary MRW 
can be obtained by just adding a friction 7 > 0, i.e., 
X At [k] = (l- 7 )X At [£:-l]-|-eAt[fc]e" At [ fe l. One can build 
a strictly increasing MRW (and consequently a stochastic 
positive multifractal measure) by just setting e^t = Ai 
in Eq. (0) and use it as a multifractal time for subordi- 
nating a monofractal process (such as an fBm). One can 
also use other laws than the (log-)normal for e and/or lu. 
Another interesting point concerns the problem of the 
existence of a limit (Ai — ► 0) stochastic process and on 
the development of a new stochastic calculus associated 
to this process. All these prospects will be addressed in 
forthcoming studies. 
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